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ABSTRACT 

We present the detections of 19 solar flares detected in high-energy gamma rays (above 100 MeV) 
"with the Fermi Large Area Telescope (LAT) during its first four years of operation. Interestingly, all 
flares are associated -with fairly fast Coronal Mass Ejections (CMEs) and are not all po"werful X-ray 
flares. We then describe the detailed temporal, spatial and spectral characteristics of the first t"wo 
long-lasting events: the 2011 March 7 flare, a moderate (M3.7) impulsive flare follo"wed by slo"wly 
varying gamma-ray emission over 13 hours, and the 2011 June 7 M2.5 flare, "which "was follo"wed by 
gamma-ray emission lasting for 2 hours. We compare the Fermi-hAT data "with X-ray and proton 
data measurements from GOES and RHESSI. We argue that a hadronic origin of the gamma rays 
is more likely than a Icptonic origin and find that the energy spectrum of the proton distribution 
softens after the 2011 March 7 flare, favoring a scenario with continuous acceleration at the flare site. 
This "work suggests that proton acceleration in solar flares is more common than previously thought, 
occurring for even modest X-ray flares, and for longer durations. 

Subject headings: Gamma rays: observations — Sun — Solar flares — Fermi Gamma-ray Space Tele- 
scope 
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1. INTRODUCTION 

Solar flares are explosive phenomena that emit elec- 
tromagnetic radiation extending from radio to 7 rays. 
It is generally agreed that magnetic energy stored in 
the solar corona and released through reconnection is 
the source of plasma heating and acceleration of elec- 
trons and ions to relativistic energies. Measurements of 
Hard X-Rays (HXRs) up to ^300 keV indicate the pres- 
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ence of electrons with energies up to few MeV producing 
bremsstrahlung in the high-density regions of the solar 
corona and chromosphere. Microwave observations indi- 
cate synchrotron emission by higher-energy, relativistic 
electrons in ~ 100 G magnetic fields. In some flares, 
often GOES [Geostationary Operational Environmental 
Satellite) X-class, electron bremss trahlung emission i s 
detected up to tens of MeV (e.g., [Trottet et al]|1998[ ). 
Nuclear 7-ray lines in the 1-10 MeV range and con- 
tinuum radiation above 100 MeV produced by acceler- 
ated protons, a particles, and heavier ions have been 
detected with instruments onboard the Solar Maximum 
Mission (SMM), the Compton Gamma Ray Observa- 
tory (CGRO)and the Reuven Rama ty High Energy Solar 
Spectroscopic Instrument (RHESSI, Lin et al.|2002 ). The 
lines are due to de-excitation of ambient (or accelerated) 
ions excited by interactions with accelerated (or ambient) 
ions. The continuum radiation is produced by interac- 
tions of >300 MeV protons and >200 MeV/n a-particles 
with ambient p rotons and Helium pr oducing neutral and 
charged pions (Murphy et al. 1987[ ). The neutral pions 
decay into a pair of 67.5 MeV 7 rays (in the rest frame 
of the pion) , and the charged pions decay ultimately into 
energetic electrons, positrons and neutrinos. The sec- 
ondary electrons and positrons emit bremsstrahlung 7 
rays in the tens of MeV energy range. These particles 
also produce inverse Compton X-rays by up-scattering 
solar optical photons and terahertz synchrotron radia- 
tion. 

In general, the 7-ray emission light curve is similar to 
that of the HXRs (possibly with some delay), lasting 
for 10 — 100 seconds. This is referred to as the "impul- 
sive" phase of the flare. However, the high-energy in- 
strument onboard CGRO Energ etic Gamma R ay Exper- 
iment Telescon e (EGRET) (Kanbach et al.||T988 , Espos^ 



ito et al. 



1999[) detected 7 rays above 100 MeV tor more 



than an "hour after the impulsive phases of 3 flares ( Ryan 
2000|. Among them, the 1991 June 11 flare is remark- 



able because the 7-ray emission (>50 MeV) lasted for 8 
hours afte r the impuls ive phase of the GOES X12.0 flare 



pCanbach et al. 1993). The measured 7-ray spectrum 



appeared to be a composite of electr on bremsstrahlung 



and pion-decay components (Kanbach et al. ||1993 Rank 
et al.||2001[ [Ramaty fc Mandzhavidzel[T99lf. The 7-ray 



light curve s h owed a. smoot h exponential decay (Kanbach 
et al. 1993). Ryan (2000) suggested that the particles 
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High-Energy Gamma-Ray Emission From Solar Flares 



accelerated during the impulsive phase of the flare could 
remain trapped for the entire duration of the flare, and 
precipitate gradually into the denser solar atmosphere to 
produce the 7 rays. Alternatively, continuous accelera- 
tion, either by a CME shock or by t urbulence in a clo sed 
magnetic loop, is a possible origin ( Rank et al.||2001 ). 

As solar activity increases with the pro gress of the so- 
lar cy cle, the Large Area Telescope (LAT, Atwood et al. 
20091) and Gamma-ray Burst Monitor (GBM, Meegan 



LAT data (Sec. 2.2 1 



et ai.l2009 ) instruments on the Fermi Gamma-Ray Space 



'1 'elescope are beginning to observe 7 rays and HXRs from 
solar flares. As we show below, in its first four years of op- 
eration, the LAT has detected emission above 100 MeV 
in at least 19 flares. In this paper, we describe the con- 
tinuous monitoring of the Sun that we perform with the 
LAT and outline broad conclusions that we derive from 
analysis of those flares. We then discuss in detail the 
first two long-duration flares: the GOES M-class flares 
SOL2011-03-07T20:12 and SOL2011-06-07T06:41. For 
both flares, the Sun was outside the field of view (FOV) 
of the LAT during the impulsive phase; nevertheless a 
significant flux of 7 rays was detected when the Sun en- 
tered the LAT FOV at intervals over the next ^ 13 hours 
on March 7 and 8 and at one interval on June 7. Results 
of the temporal and spectral analyses and localization 
studies are described in Sec. [3] followed by a brief inter- 
pretation in Sec. |4] 

2. DATA ANALYSIS 

The Fermi-LAT is a wide field of view (FOV), imag- 
ing telescope for high-energy 7 rays, designed to cove r an 
energy range from 20 MeV up to more than 300 GeV ( At- 
wood et al |2009|) . The instrument consists of a precision 



ige tro 
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tracker with silicon strip detectors above a cesium-iodide 
calorimeter. Both are enclosed in the plastic scintillators 
of the Anti- Coincidence Detector (ACD) that provides 
charged-particle tagging for background rejection. For 
bright solar flares an intense flux of X-rays during the 
impulsive phase of the flare can result in pulse pile-up 
in the ACD scintillators within the integration time of 
the ACD readout. A coincident 7 ray entering the LAT 
within that integration time can be misidentified by the 
instrument's flight software or event-classification ground 
software as a charged particle and thereby mistakenly ve- 
toed. During these periods, the nominal LAT instrument 
response functions do not apply, and the data cannot be 
analyzed by stan dard software. These issues were ad- 



dressed in detail in Ackermann et al. (2012b) for the 2010 
June 12 flare. The LAI' instrument team closely moni- 
tors for this effect and tags such data as "bad" through 
redeliveries to the public data archivcp°| Data marked 
as "bad" have been included in the continuous monitor 
of the Sun that we describe in Sec. |2.1[ The goal of this 
monitoring is to detect any possible increase of the solar 
flux. (The effect of the ACD pile-up is to decrease the 
effective collecting area; therefore a significant excess will 
still correspond to an increase of the flux) . On the other 
hand the instrument response functions (used for high- 
level analysis) does not account for this effect, and the 
value of the measured flux will not be correct during time 
intervals with high ACD pile-up. We therefore remove 
such time intervals in the detailed analysis of standard 



2.1. LAT SunMonitor 

Ferm,i has spent more than 95% of its mission to date 
in survey mode, in which the spacecraft rocks to put the 
center of the LAT FOV 50° north and 50° south of the 
orbital equator on alternate orbits. In this way the LAT 
monitors the entire sky every two orbits, or about every 
three hours, and observes the Sun for ~20 — 40 con- 
tiguous minutes in that time (see Sec. 3.1). We have 



created an automated data analysis pipeline, the Fermi- 
LAT SunMonitor, to monitor the high-energy 7-ray flux 
from the Sun throughout the Fermi mission. The time 
intervals during which we run the analysis are the in- 
tervals in which the Sun is less than 60° off-axis for the 
LAT. In this way, each interval corresponds to the max- 
imum time with continuous Sun exposure, and the du- 
rations of these intervals vary as the Sun advances along 
the ecliptic and as the orbit of Fermi precesses. We use 7 
rays with energies betwe en 100 MeV and 10 GeV f rom the 
P7S0URCE_V6 data class ( [Ackermann et al.|2012a| , which 
is well-suited for point-source analysis. Contarnination 
from 7 rays produced by cosmic-ray interactions with 
the Earth's atmosphere is reduced by selecting events 
measured to be within 105° of the zenith. Each inter- 
val is analyzed using a Region Of Interest (ROI) of 12° 
radius, centered on the position of the Sun at the cen- 
tral time of the interval. The maximum deviation of the 
true position of the Sun during these ^ 30 minutes due 
to its apparent motion is less than 1', which is smaller 
than the typical angular resolution of the instrument (the 
68% containment angle of the reconstructed incoming 7- 
ray direction for normal incidence at 1 GeV is 0?8 and 
at 100 MeV is 6°) and than the localization precision for 
even bright solar flares. It is therefore not necessary to 
apply a correction to account for the motion of the Sun 
from the center of the ROI. 

In each time window, we perform a spectral analy- 
sis using the unbinned maximum likelihood algorithm 
gtlika^ The ROI is modeled with a solar component 
and two templates for diffuse gamma-ray background 
emission: a Galactic component produced by the inter- 
action of cosmic rays with the gas and interstellar radia- 
tion fields of the Milky Way, and an isotropic component 
that includes both the contribution of the extragalactic 
diffuse emission and the residual cosmic rayfQ We fix 
the normalization of the Galactic componentbut leave 
the normalization of the isotropic background as a free 
parameter. 

We verified that the background model describes the 
data well at times away from flares. To test for tran- 
sient solar emission, the Sun is assumed to be a point 
source with a 7-ray spectrum described by a power law 
with an exponential cut-off (see Sec. 3.3). Th e three 

" " we per- 



3.3 



^^ http : //f ermi . gsf c . nasa . gov/ssc/data/ access/ 



spectral parameters are left free. (In Sec. 
form spectral analyses for two flares using an additional, 
physically motivated spectral form.) In fitting all free 
parameters to maximize the likelihood, we compute the 

^^We used the ScienceTools version 09-28-00, available on the 
Fermi Science Support Center web site http://fermi.gsf c.naseiT] 
gov/ssc/) 

"^'i'he models used for this analysis, gal_2y6arp 7v6-v0.f Its and 
iso_p7v6source.txt, arc available at http://f erml.gsf c .nasa. | 
Igov/ssc/data/access/lat/BackgroundModels.htmTl 
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significance of the point source using the Test Statistic, 
TS=2[log£ — log£o], where £o is the likelihood of the 
null hypothesis (no source present at the position of the 
Sun) and C is the maximum likelihood when the source 
is added to the model. 

To understand the statistical significance correspond- 
ing to a particular value of TS, we applied the 
SunMonitor analysis to a test location moving along 
the ecliptic plane 180° from the Sun for the full four- 
year data set (corresponding roughly to 10^ realizations). 
Over such a long period, this fiducial location samples 
the same charged-particle and celestial backgrounds as 
the Sun but is free of any possible flare signal. The dis- 
tribution of TS values determined for the test location 
was consistent with 1/2x2- From Wilks's theorem, we 
might naively expect, with the addition of three param- 
eters (flux normalization, photon index, and cutoff en- 
ergy), that TS would be distributed as xi- In our case, 
however, the factor of ^ arises from the requirement that 
source flux not be negative, and the reduction in degrees 
of freedom from 3 to 2 results from correlation among 
model parameters. 

Continuous monitoring of the Sun has led to the high- 
confidence detections of a number of flares with the LAT. 
In Table [TJ we list detections for which the TS is greater 
than 30 (roughly corresponding to 5 cr), along with es- 
timates of their >100MeV average flux during the in- 
dicated durations. We have grouped consecutive detec- 
tions together into 16 distinct flaring episodes. To com- 
plement the results of the Sun Monitor, we also ana lyzed 
the LAT Low Energy (LLE) ( [Pelassa et "aLl|2010[ ) data 
for every flare detected by Fermt-iJiiM, detecting three 
additional flares. These data have relaxed event selection 
compared to the data class used by the SunMonitor and 
are not compromised by pileup effects in the ACD dur- 
ing the impulsive phase of a flare. LLE analysis achieves 
larger effective area to transients in the 30 MeV to 1 
GeV band than standard analysis. We have indicated 
the LLE detections during the impulsive phase using the 
label "LLE" in the TS column. Although it is pos sible to 
measure fluxes usin g LLE data (see, for example, Acker- 
mann et al.|2012b" ), a dedicated spectral analysis of each 
flare would be necessary. This is beyond the scope of this 
paper, and we defer to further publications for evaluation 
of the fluxes of LLE flares. CME shock spe eds measured 
by SOHO-LASCO ( |Brueckner e"raL]|l995] ) are obtained 
from the SOHO/LASCO CME online catalog and are 
reported in the third column of the table. 

Several features are immediately apparent from the ta- 
ble. Although the SunMonitor has analyzed data since 
the start of the mission, the detection of solar flares starts 
only in mid 2010, with the rise of solar activity in the 
current cycle. In flve cases, labeled as type "I" for "Im- 
pulsive" in the table, we detected the impulsive phase 
above 30 MeV using LLE data, and in two of them the 
flare was also detected above 100 MeV in standard like- 
lihood analysis by the LAT SunMonitor. In some cases, 
due to the partial overlap between the duration of the 
GOES X-ray pulse and the time interval in which we 

^^The CME catalog is generated and maintained at the CDAW 
Data Center by NASA and The Cathohc University of America 
in cooperation with the Naval Research Laboratory. SOHO is a 
project of international cooperation between ESA and NASA. 



detect the solar flare at high energy, we cannot disen- 
tangle the impulsive phase from the long-duration (or 
"Sustained") emission, and we have labeled these cases 
as type "I/S". When high-energy emission (>100 MeV) 
is detected in time intervals subsequent to the impulsive 
HXR emission, we have labeled the detections as type 
"S" for long-duration, "Sustained" emission. In three 
cases (2012 January 27; 2012 March 09; 2012 March 10) 
the Sun was in the LAT FOV at the time of the HXR 
impulsive flare, but the LAT did not detect it, suggesting 
that high-energy emission can also arise at later times. 
We have labeled these type "D" for "Delayed." Although 
these cases are rare, they are particularly interesting and 
will be the subject of further analysis in a subsequent 
paper. 

In almost all the cases, the flares are associated with 
moderately bright X-ray flares, although unlike the 
EGRET-detected flares, not all are X-class. Instead, 
they are predominantly M-class and, in one case (2011 
June 02), high energy emission occurred coincidently 
with a series of C-class flares. All the flares are Solar 
Eruptive Events (SEE), i.e., they are associated with 
CMEs and SEPs. Most have fast (>500 km s^^) CMEs, 
and six have CME velocities ^ 2000 km s^^. The total 
energy radiated in > 100 MeV 7 rays varies from < 10^^ 
to 10^^ ergs and thus spans a broader dynamic range 
than other characteristics (e.g. GOES flux spans two 
decades, CME speed a factor of 4, etc.). It is nonethe- 
less small compared to the typical total energy of a flare, 
which can be > 10^^ ergs for the largest flares, attesting 
to the high sensitivity of the Fermi-LAT . We caution 
that flares whose only high-energy emission is during the 
impulsive phase can be missed because of the modest 
in-aperture viewing fraction afforded by the sky survey 
observing strategy. Averaged over year, the Sun is within 
the LAT FOV --20% of the time, and only a small num- 
ber of flares is viewed during their impulsive phase. An 
example of this is the X5.4 flare on 2012 March 7 whose 
main impulsive burst was missed by Fermi. 

2.2. Detailed Analysis of Two Flares 

A detailed and complete analysis of all t he flares in 
Table fll is be yond the scope of this paper. 
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et al. ([2012b ) discuss a flare that exhibited only impulsive 
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emission (<1 minute in duration) at high energies, the 
M2.0 flare of 2010 June 12, which showed 7-ray lines 
below 10 MeV and continuum up to 300 MeV in the LAT. 
Here we focus on the first two long-duration fiares: the 
M3.7 solar flare of 2011 March 7, which wa s detected with 
LAT over a interval of almost 14 hours ("Allafort et al 
2011| , and the M2.5 flare of 2011 June 7 (Ta naka eTal 
2011 1 detected for less than an hour. For these two flares 
we used P7S0URCE_V6 event class data and instrument 
response functions in the 60 MeV to 6 GeV energy range 
for spectral analysis and the 100 MeV to 6 GeV range 
for localization. We restricted the data set to 7-rays 
arriving with zenith angles less than 100° to minimize 
contamination from atmospheric 7 rays, and we analyzed 
7-rays within a 12° ROI around the Sun. We included the 
azimuthal (0) dependence of the effective area of the LAT 
when calculating the exposure for the likelihood analysis. 
While the (f) dependence averages out for observing time 
scales of days and longer, on scales of minutes and hours 
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TABLE 1 
Solar flares detected by Fermi LAT from 2008 August to 2012 August. 



Date (UT) 


Duration 
mill. 


GOES X-ray 
Class, Start-End 


OMEt 
Speed, km s^^ 


TS 


Type 


Flux (>100 MeV) 
xlO^^ph cm-2 s-i 




2010-06-12 00:55 


~1 


M2.0, 


00:30-01:02 


486 


LLE* 


I 


(-) 


2011-03-07 20:15 
23:26 

2011-03-08 02:38 
05:49 


25 
36 
35 
35 


M3.7, 


19:43-20:58 


2125 


230 
520 
450 
200 


1/S 
S 
S 
S 


(1.9± 0.3) 
(3.5± 0.3) 
(3.5± 0.3) 
(1.9± 0.3) 


2011-06-02 09:43 


45 


02.7,9:42-9:50 


976 


35 


1/S 


(0.4±0.2) 


2011-06-07 07:47 


53 


M2.5, 


06:16-06:59 


1255 


570 


S 


(3.6 ± 0.3) 


2011-08-04 04:59 


34 


M9.3, 


03:41-04:04 


1315 


390 


S 


(2.5 ± 0.3) 


2011-08-09 08:01 


^ 1 


X6.9, 


07:48-08:08 


1610 


LLE* 


I 


(-) 


2011-09-06 22:17 
2011-09-06 22:13 


< 1 
35 


X2.1, 


22:12-22:24 


575 


LLE* 
2600 


I 

I/S 


(-) 
t 


2011-09-07 23:36 


63 


X1.8, 


22:32-22:44 


792 


350 


s 


(1.0 ± 0.1) 


2011-09-24 09:35 


~1 


X1.9, 


09:21-09:48 


1936 


LLE* 


I 


(-) 


2012-01-23 04:07 
05:25 
07:26 
08:47 


51 
69 
16 
35 


M8.7, 


03:38-04:34 


1953 


180 
650 
69 
97 


1/S 
S 
S 
S 


(0.8 ± 0.1) 
(2.1±0.2) 
(3.7±0.9) 
(2.6±0.5) 


2012-01-27 19:45 
21:13 


11 

24 


X1.7, 


17:37-18:56 


1930 


78 
47 


D 

S 


(3.2±0.8) 
(l.OitO.3) 


2012-03-05 04:12 
05:26 
07:23 


49 
71 

28 


XI. 1, 


02:30-04:43 


1602 


69 
250 
39 


1/S 
S 
S 


(0.5±0.1) 
(0.9±0.1) 
(0.8±0.2) 


2012-03-07 00:46 

2012-03-07 03:56 
07:07 
10:18 
13:29 
19:51 


31 

32 
32 
32 
32 

25 


X5.4. 
X1.3, 


00:02-00:40 
01:05-01:23 


1785 


22000 

16000 

8900 

1900 

120 

50 


S 

I/S 
S 
S 
S 
S 
S 


t 

(113.1±2.0) 
(71.9±1.6) 
(30.1±1.5) 
(8.9±1.9) 
(0.4±0.1) 


2012-03-09 05:17 
06:52 
08:28 


34 
35 

34 


M6.3, 


03:22-04:18 


844 


51 
100 
159 


D 

S 
S 


(0.6±0.2) 
(0.9±0.2) 
(1.4±0.2) 


2012-03-10 21:05 


30 


M8.4, 


17:15-18:30 


1379 


43 


D 


(0.4±0.1) 


2012-05-17 02:18 


22 


M5.1, 


01:25-02:14 


1582 


45 


1/S 


(l.OitO.3) 


2012-06-03 17:52:33 

17:40 


~1 
23 


M3.3, 


17:48-17:57 


605 


LLE* 
300 


I 

I/S 


(-) 
(3.2±0.4) 


2012-06-14 14:48 


49 


M1.9, 12:52-15:56 


987 


49 


I/S 


(1.1±0.3) 


2012-07-06 23:19 


52 


XI. 1,23:15-23:49 


892 


930 


I/S 


(3.5±0.2) 


T CME data are available at the following 
^ The flux estimate is unreliable because 
* LLE detections arc >30 MoV while TS 


url: http;//cdaM. gsfc.nasa.gov/CME_list/ 
of X-ray pile-up in the ACD. 
values are calculated for >100 MeV. 









- those relevant to solar flare analysis - the range of (f> 
angle for an individual source is not well-represented by 
the azimuthal average: e.g., the effective area can differ 
from the azimuthal averaged by 5% typically and more 
than 10% below 100 MeV or far off-axis (incidence angle 

>6o°g 

It is important to account for the apparent motion of 
the Sun in the analyses of particularly long flares, i.e. 
when the motion is a non-negligible fraction of the lo- 
calization accuracy which is a tighter constrain than the 
size of the LAT point spread function (PSF). Such is 
the case in analyzing the ~14 hours of 7-ray emission 
from the March 7 event as a whole. (It is not necessary 
for the June 7 event since the high energy emission was 
visible by LAT for only 36 minutes.) We have devel- 
oped a dedicated "Sun-centering" analysis tool for mov- 
ing sources. We transform the directions of all the 7 rays 
into ecliptic coordinates, then translate them in ecliptic 
longitude to keep the Sun at position (0,0) as time passes. 



® http : //f ermi . gsf c . nasa . gov/ssc/data/analysis/ 
ldocumentation/Cicerone/Cicerone_LAT_IRFs tor more details. 



We apply the same operation to the pointing history of 
the LAT to keep an accurate account of the exposure. 
In Sun-centered coordinates the diffuse backgrounds are 
well approximated by an isotropic intensity. We verified 
that fitting the standard isotropic template to the data 
with a free normalization coefficient is provides a good 
representation of the backgrounds. Two additional back- 
grounds releva nt for analyses of s olar flares are the quiet 
Sun emissions ( Abdo et al.||20lf ) from cosmic-ray inter- 
actions with the Sun (disk) and with the solar radiation 
field. For the typical duration of the SunMonitor analysis 
(nominally less than an hour), these backgrounds are not 
significant. For long duration flares (such as the March 
7), or for detailed localization analysis, it is important to 
account for the solar disk component, modeled as a point 
source at the position of the center of the d isk with the 
param eters flxed to their measured values (Abdo et al. 



2011|. 

Uncertainties in the calibration of the LAT introduce 
systematic errors on the measurements. Effective area 
uncertainty is dominant, and for the P7SDURCE_V6 event 
class it is estimated to be ^10% at 100 MeV, decreasing 
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to -5% at 560 MeV, and increasing to ~10% at 10 GeV 
and above. In order to estimate the systematic uncer- 
tainties in the model parameters, we repeat the analysis 
using a set of custom modified instrument response func- 
ti ons. We describe this "bra cketing" technique in detail 
in 



Ackermann et al. (2012a) 



3. RESULTS 

On 2011 March 7, the solar activity increased dramat- 
ically, with a dozen M-class solar flares detected during 
the subsequent two days with the GOES soft X-ray mon- 
itor. Intense HXR emission (up to 300 keV) observed by 
RHESSI accompanied an M3.7 solar flare that erupted 
from the NOAA Active Region (AR) 11164 in the north- 
west quadrant. The impulsive phase started at about 
19:43 UT and ended at 20:10 UTHAround 6:16 UT on 
2011 June 7, an M2.5 flare erupted from AR 11226 in the 
southwest quadrant, with prominent soft X-ray and HXR 
emission, ending at 06:59 UTpjThe heliographic longi- 
tude of the flare site was similar to that of the March 7 
flare. For both flares, the impulsive phase seen in HXRs 
by RHESSI occurred entirely while the Sun was outside 
the Fermi-LAT FOV. Without any interference from the 
impulsive X-rays in the ACD the Fermi-LAT was able 
to start observing the Sun 33 and 92 minutes after the 
start of the March 7 and June 7 events, respectively. 

3.1. Gamma-ray Light Curve 

We used the output of the SunMonitor to identify the 
intervals during which significant emission was detected 
for the March 7 and June 7 flares. We then performed 
a more detailed spectral analysis for these intervals and 
the preceding and following six-hour intervals. The up- 
per panels of Figures [l] and [2] show the flux measure- 
ments or upper limits in each observing window, where 
each window is defined to be when the 12°-radius anal- 
ysis ROI (see Sec. [2|around the Sun is entirely within 
the 70° FOV. Table[2] gives the precise time and dura- 
tion of each window in which we found a positive de- 
tection of emission at energies greater than 100 MeV. 
(Since the flare is detected with high significance in four 
consecutive time windows, we include also the last point 
with TS=20.) The most striking feature of the March 
7 flare is the slow increase of the flux, which reaches its 
peak, F{> lOOMeV) = (3.8 ± 0.3) x 10"^ ph s'^ cm'^, 
more than 7 hours after the onset of the M3.7 flare it- 
self. The flux then decreased gradually until the last 
significant detection, almost 14 hours after the impul- 
sive phase. The June 7 flare was detected significantly 
only in the first observing window after the flare on- 
set, with only upper limits obtained for the following 
windows. The fiux measured in this first window, less 
than two hours after the impulsive phase, is of the same 
magnitude as the peak flux for the March 7 flare, with 
F(> lOOMeV) = (3.4 ±0.2) x 10"^ ph cm-^ s"!. The 
upper limits shown before and after the detections are 
consistent with the quiet Sun flux of (4.6 zt 0.2) x 10~^ 
ph cm~^ s~^ reported in Abdo et al. (2011 ). 

We also searched for spectral evolution from one ob- 
serving window to t he next for the March 7 flare, as 
described in Sec. I 



3.2. Localization 

Localization of the high-energy 7-ray emission relative 
to the X-ray flaring site can provide valuable information 
about the source of the accelerated particles producing 
the high-energy emission. The key question is whether 
that emission is tightly constrained to the flaring site, 
or displaced, or broadly distributed. This is complicated 
by the fact that on March 7 and 8, within the 10 hours 
following the M3.7 flare from AR 11164, other regions of 
the Sun were active (AR 11165 and 11171), producing a 
total of three M-class flares (Figure [I]). It is necessary 
to know the location of the source or sources of the high 
energy emission to interpret the 7-ray light curve prop- 
erly. We used the maximum likelihood fitting package, 
pointlike, which was used in the localiza tions for the 
sourc es in the LAT 2-year catalog (2FGL, Nolan et al. 
2012 1 sources and is particularly suited for simultane- 
ously fitting position, spec truin and possible spati al ex- 
tension of a source (Kerr||2011 Lande et al. |2012 ). We 
limited this analysis to energies greater than 100 IVleV to 
exclude 7 rays with the most uncertain directions. 

For both flares we modeled the long-duration emission 
with a point source and found the best-fit source loca- 
tion for the full duration of the high-energy detection. 
The best locations and corresponding circular 68% and 
95%-confidence statistical uncertainty regions are plot- 
ted for both flares on images of the Sun from the Atmo- 
spheric Imaging Assembly (AIA) onboard the Solar Dy- 
namic Observatory (SDO) (Figure l3|. Positions are ex- 
pressed in Cartesian-projected heliographic coordinates. 
For the March 7 event, the best position is (750" , 690" ) 
offset from the center of the solar disk with 68% and 
95%-confldence error radii of 290", and 500", respec- 
tively. This position in the northwest quadrant is consis- 
tent with that of AR 11164 as imaged by RHESSI and 
SDO. For June 7, the best location is (870", -350") with 
a 68%-confidence radius of 500", in the southwest quad- 
rant of the solar disk and consistent with the position of 
AR 11226. The uncertainty is large because the source 
flux is modest. 

A possible bias in the determination of the position 
of an off-axis 7-ray source, particularly at low energy, is 
the so-called "flsheye" effect, a known systematic shift in 
the reconstruction of each event toward the center of the 



LAT FOV (Ackermann et al.||2012a). This bias does not 



http : //www . swpc . noaa . gov/weekly/pdf /prf 1854 . pdf I 
^^Ettp : //www . swpc . noaa . gov/weekly/pdf /prf 1867 . pdf 



affect long observations of steady sources with the LAT, 
but for short observations at high incidence angle, as is 
the case for most solar flares, this effect can be impo rtant, 
e.g., as we reported f or the 2010 June 12 solar flare ( Ack- 
ermann et al.|2012b| ). From Monte Carlo simulations we 
evaluated the energy-dependent correction to be applied 
to the incidence angle of each 7 ray, and for the March 7 
flare the average shift of the reconstructed direction with 
respect to the Monte Carlo position is ~100". This bias 
is negligible compared to statistical uncertainties, so we 
have not included it in the best-flt positions we reported 
above. 

For the March 7 flare, we repeated the localization 
analysis for each observing window with a positive de- 
tection, but the low counting statistics produced large 
uncertainties that are not constraining, and much larger 
than the solar disk. We were therefore unable to evaluate 
whether the high-energy emission site moved during the 
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Fig. 1.— Multi-wavelength and proton light curves of GOES M3.7 SOL2011-03-07T20:12. Gamma-ray data are from the Fermi-LAT, 
HXRs from RHESSI, soft X-rays and protons from GOES. Vertical error bars of LAT data indicate la statistical uncertainties, and gray 
arrows are 95% upper limits. The horizontal bar for each flux point represents the true duration over which the flux was computed (see 
Table [2]|, the Sun being out of the FOV at other times. 



event. 

Using pointlike we investigated the possibility that 
the 7-ray emission is spatially extended and found no 
significant TS increase for an extended source relative to 
a point source, for either the total time of detection or 
for individual intervals. 

3.3. Spectrum 

We can evaluate the contributions from two distinct 
emission mechanisms - bremsstrahlung from accelerated 
electrons, and decay of pions from interactions of accel- 



erated hadrons - regardless of how the charged parti- 
cles are accelerated by comparing the > 100 MeV data 
to model predictions. We used a simple power law 
dN/dE = NqE~^ , where F is the photon index and E 
the energy of the 7-rays, to describe bremsstrahlung from 
a non-thermal electron distribution. To test whether the 
electron distribution breaks or cuts off at high energies, 
which would be reflected as a break or cut-off in the 7- 
ray distribution at lower energies, we also considered a 
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Fig. 2.— Multi-wavelength and proton light curves of GOES M2.5 SOL2011-06-07T06:41. Gamma-ray data are from the Fermi-LAT, 
HXRs from RHESSI, soft X-rays and protons from GOES. Vertical error bars of LAT data indicate la statistical uncertainties, and gray 
arrowrs are 95% upper limits. The horizontal bar for each flux point represents the true duration over which the flux was computed (see 
Table [2]|, the Sun being out of the FOV at other times. 



power law with an exponential cut-off (ExpCutoff): 

dN{E) 



dE 



NoE- 



(1) 



Although the two models are identical only for E^o —> 
oo and the likelihood-ratio test cannot be rigorously ap- 
plied, we report in Table [3] the increment in TS obtained 
by including the exponential cut-off in the model. We 
also report the best-fit values of the parameters of the 



two models for both flares. 

The power law with an exponential cut-off is a better 
representation of the 7-ray spectrum as the increment 
in TS value is greater than 25 in both cases. This sim- 
ple analytic function is very similar to the 7-ray spec- 
trum resulting from the decay of pions. For the pion- 
decay model, we used calculated 7-ray spectra resulting 
from interactions of protons and alpha particles having 
an isotropic momentum distribution and a power-law ki- 
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TABLE 2 
Fermj-LAT Observing windows, duration, 7-RAY flux, and best-fit 

PROTON SPECTRAL INDEX. 



Date (UT) 



Duration Flux (>100 MeV) 
min. xlO~^ ph cm~^ s^'^ 



Proton index 



GOES M3.7 flare, SOL2011-03-07T20:12 



2011-03-07 20:15:42.6 
2011-03-07 23:26:51.6 
2011-03-08 02:37:37.6 
2011-03-08 05:49:03.6 
2011-03-08 09:13:06.7 



24 


1.7±0.2I^;^ 


4.0±0.5I^;^ 


33.5 


3.3±0.3l°-^ 


4.6±0.3tol 


34 


3.5±0.3+°:^ 


4.9±0.3t°-^2 


34 


1.8±0.2+0-2 


>5.6 


21 


0.4±0.liro'3 


t 



GOES M2.5 flare, SOL201 1-06-07X06:41 



2011-06-07 07:47:40 



36 



3.1±0.2lj!;j 



4.3±0.3: 



In this time interval, tlic number of -y-rays is small and the pion-deeay template spcetrt 
.oes not produee a statiatieally satisf'aetory fit. The best-fit model to the "y-ray data 
escribed by a power-law with spectral index r = 2.7±0.4 with the reported flux. 



SOL2011-03-07T20:12 

O. Ra.Dec^ 347.031 -5.557 l^.B^.Pr 119.4 -7.2 -22.8 



SOL2011-06-07T06:41 

Q: Ra,Dec; 74.908 22.712 l^.B^.P: 159.8 0.0 -13.1 



1000 



500 



y 



-500 



-1000 




1000 



500 



E 



-500 



-1000 




-1000 



-500 



X [arcsec] 




X [arcsec] 



500 



1000 



Fig. 3. — Fermi-LAT localization of the high-energy emission (100 MeV-10 GeV) for the 2011 March 7 (left) and June 7 (right) flares. 
In each image, the Sun is displayed in heliographic coordinates, where the projection of the rotational axis of the Sun is along the Y-axis, 
the Z-axis is the line of sight (from the Sun to the observer) and the X-axis is the remaining orthogonal basis vector. The grid of constant 
heliographic latitude and longitude is yellow and lines of constant R.A. and Dec. (J2000) are white. The backgrounds are SDO/AIA images 
at 304 A of the chromosphere and transition region around the Sun at the given epoch. The labels indicate the NOAA numbers of the 
active regions. The LAT localizations of the flares are shown by the green circles, solid and dashed lines at the 68% and 95% error radii 
respectively. 



netie energy s pectrum w i th ind ex s 

(11987) 



Np{£) oc e *] as de- 
We obtained a series of 



scribed in Murphy et al 
tabulated photon spectra by varying the index s of the 
protons. We compared the associated pion-decay photon 
spectra with the data by performing an unbinned likeli- 
hood spectral analysis and computing the value of the 
logarithm of the likelihood [log£] with the normaliza- 
tion of the photon spectral template as a free parameter. 
To estimate the maximum likelihood proton index and 
its statistical uncertainty, we fit the values of — log C cal- 
culated for the tabulated model spectra near their min- 
ima with a parabolic function of proton index s. The 
minimum gives the most likely index s = sq for the pion- 



decay model. We calculate the statistical uncertainty 
on the proton index using ATS= 2[logCs — log£o]=li 
which corresponds, by Wilks's theorem, to the 68% Con- 
fidence Level (CL) for a x^ distribution with one degree 
of freedom. In the last row of Table [Sj we report the 
best-fit values of the proton spectral indices and their es- 
timated uncertainties for both flares. Systematic errors 
are estimated using the b racketing method described in 



Ackermann et al. 



J2012a |. 
I'igure |3] compares the observed and predicted num- 
ber of 7-ray counts for different models for the March 
7 and June 7 solar flares, while Figure [5] shows the 
reconstructed spectral energy distribution for the two 
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TABLE 3 

Best-fit parameters with statistical and systematic errors for the 

TIME integrated SPECTRA. 



Model 






SOL2011-03-07T20:12 


SOL2011-06-07T06:41 


Power Law 




Fluxt 

r 


L88±0.05±0.2 
2.56±0.06±0.08 


2.6±0.2±0.2 
2.45±0.09±0.08 


ExpCutoff 


Eco 


Fluxt 

r 

[MeV] 
ATS* 


2.1±0.1±0.2 

1.5±0.3+_lil 

130±20li''° 

96 


3.0±0.2+^-^ 
1 1+0 4+0.1 

210±40l°o* 
44 


Pion decay 


Fluxt 
Proton index 


2.1±0.1±0.2 
4.5±0.2±0.2 


3.1±0.2+j;-j 
4.3±0.3l|]J 



The integral flux bet 



LOO MeV and 10 GeV is in 
nt for the exponential cut-off model : 



flares. To obtain the 10 model-independent photon spec- 
tral data points, we fit the background through the en- 
tire energy range (60 MeV to 6 GeV), then within each 
energy bin, determined the source flux using the Fermi 
Science Tool gtlike assuming a power-law photon spec- 
trum having fixed index of 2 and background parameters 
fixed at those obtained from the fit to the entire energy 
range. We note that the likelihood spectral fltting does 
not account for the finite ene rgy resolution of the LAT 



15% at 100 MeV; Ackermann et al. 2012a 



Atwood et al. 2009). The uncertainties on the energy 



measurement are smaller than the width of the spectral 
feature (curvature) we measure and therefore do not sig- 
nificantly affect our results. 

Figure [S] shows the detections (TS>9) and 95% CL 
upper limits that result from this analysis. Statistically 
significant solar emission up to ~1 GeV is apparent for 
the March 7 fiare. Also shown in Figure [5] are the best- 
fitting photon spectral forms for the pion-decay, expo- 
nential cut-off, and power-law models. The insets show 
the likelihood profiles of the pion-decay model as a func- 
tion of the proton spectral index. 

For the March 7 flare, the best-flt index for a power- 
law photon model is F ^ 2.56±0.06, but this clearly 
does not provide a good flt to the data. Including 
the exponential cut-off makes a statistically significant 
improvement, with a ATS — 96. In this case, the 
power-law component has an index of F ^ 1.5±0.3 and 
Eco ~ 130±20 MeV. It is apparent from Figure [5] that 
the exponential cut-off model with these parameters is 
similar in shape to the pion-decay model. The best fit 
for the March 7 flare is obtained with the pion-decay 
model with proton index s ^ 4.5±0.2. Integrating the 
best-fitting photon spectrum from 100 MeV to 10 GeV, 
the flare-averaged 7-ray flux is F(lOOMeV-lOGeV) -- 
(2.1±0.1)xl0~^ ph cm~^ s~^, which corresponds to an 
energy flux of (7.2±0.4)x 10~^ ergcm"^ s~^. For the du- 
ration of the 7-ray emission, about 13.3 hours, the total 
emitted energy > 100 MeV is £^ ^ (7.7±0.4)x 10^^ erg. 

For the June 7 flare, the best-fit photon index for the 
power-law model is F ^ 2.45±0.09, but the exponen- 
tial cut-off model is preferred (with a ATS = 44). In 
this case, Eco ^ 210±40 MeV. For the pion-decay model, 
the best-fit proton index s ~ 4.3±0.3 suggests a slightly 
harder population of accelerated protons compared to the 
March 7 flare. Using this model, the average 100 MeV to 



TABLE 4 

Spectral energy distributions for the 

March 7 and June 7 solar flares. 



Energy Bin 
MeV 



60-95 

95-150 

150-239 

239-378 

378-600 

600-950 

952-1508 

1509-2391 

2391-3789 

3780-6000 



Flux 

xlO~^ergcm~^ s~^ 

2011 March 7 2011 June 7 



2.7 ± 0.4 

3.8 ± 0.4 
4.5 ± 0.4 
4.2 ± 0.4 
1.7 ± 0.3 
0.6 ± 0.2 

< 0.5 

< 0.7 

< 0.6 

< 0.9 



2.5 ± 0.8 
5.3 ± 0.9 
7.7 ± 1.0 

4.6 ± 0.8 

3.7 ± 0.9 

< 2.1 

< 1.3 

< 1.9 

< 4.1 

< 4.0 



10 GeV flux is (3.1±0.2)xl0"''^ phcm^^ g-i^ correspond- 
ing to an energy flux of (10±l)xl0~^ ergcm"^ s~^. The 
total emitted energy > 100 MeV over the 36 minutes of 
detection of the flare was £^ ~(4.8±0.5)xl0^-^ erg. 

The March 7 event was sufficiently strong that we could 
study spectral evolution. We carried out the spectral 
analysis method for each of the five time intervals in Ta- 
ble[2] For no time interval other than the last one, where 
the emission is weak, was the simple power law model an 
adequate fit to the data. The first four time intervals 
are well described by the pion-decay model and show 
a significant, monotonic softening in time of the proton 
spectrum responsible for the 7-ray emission. The best-fit 
proton spectral indices and the corresponding integrated 
100 MeV to 10 GeV fluxes are given in Table [2] 

3.4. Multi-wavelength and Proton Data 

Both fiares were observed with instruments in X- 
rays. Extreme Ultraviolet (EUV), and radio as well as 
by charged-particle detectors. They were both modest 
GOES M-class flares with no RHESSI HXR signal above 
300 ke\pj and both were associated with SEE. The light 
curves oTthe GOES 0. 1-0.8 nm band, the RHESSI 25-50 
keV channel, and GOES proton data are plotted below 
the Fermi-LAT light curves in Figures [l] and [2] For the 
March 7 event, we note additional peaks in theGOES X- 
ray light curves that correspond to flares from ARs other 



Hi 



^http://sprg.ssl.berkeley.6du/-tohban/browser/?show= | 
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Fig. 4. — Comparison between observed 7-ray counts and predicted 7-ray counts for different models for the March 7 and June 7 solar 
flares. The black solid line is the histogram of the observed counts. Vertical bars are plotted using the Gehrels (1986) formula to show 
the expected uncertainty associated to each bin. Red, blue, and green solid lines are the expected numbers ot counts tor different models 
(power-law, power-law with expone ntial cut-off, and pion-decay, respectively) after likelihood maximization. The contribution of the quiet 
Sun is shown as green dashed line ||Abdo et al.|201l)|. 



■K° ,7r^ decay 
Power law 
Cut-off Power law 






— 7r",7r^ decay 

- - Power law 




j Cut-off Power law 




^'"^^ 1 


■ — 


- 




+ X 


~ "^ - ^ — 


,— 












t- 








3.0 


-2Alos(£) 




\ 
















2.5 




+ \ 






2.0 






\ 


1.0 


I ,-' 


\ 


0.5 

0.0 

4 




\ 


1,2 AA .1,6 4,8 5 
Proton Spectral Index 









\ 







10^ 
Energy [MeV] 



10" 
Energy [MeV] 



Fig. 5. — Reconstructed Fermi-hAT spectral energy distributions of tlie March 7 and June 7 solar flares. Tlie red solid, dotted, and dashed 
lines represent the pion-decay, exponential cut-off, and power-law models, respectively. Vertical bars indicate la statistical uncertainties. 
Upper limits are computed at the 95% CL. The insets show the profile of the statistic ATS=— 2A log £ as a function of the proton spectral 
index used to calculate the pion decay spectrum. The horizontal dashed lines at 1 indicate the value used to compute the uncertainty on 
the proton index corresponding to the 68% CL for a x^ distribution with 1 degree of freedom. Parameter values for the models obtained 
by maximizing the likelihood and used to evaluate the spectral fits shown above are in Table |3] flux measurements with statistical errors 
are reported in Table |4] 
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Fermi-ljAT collaboration 



than that producing the M3.7 flare: AR 11166, 11165, 
11171. For the M3.7 flare starting at 19:43 UT on March 
7, the impulsive phase lasted about 15 minutes, and the 
HXR and soft X-ray emissions have light curves and 
intensities typical of many flares. The bremsstrahlung 
HXR emission indicates acceleration of electrons at least 
to 1 MeV. The total energy in 20-300 keV HXRs inte- 
grated over the impulsive phase is Shxr ^1.7x10^^ erg, 
--200 times larger than the energy £j{E > lOOMeV) 
inte grate d over the 13.3 hours of 7-ray emission (from 
Sec. 3.3 1. RHESSI imaging of the March 7 flare shows 
an unusually long (^ 10^° cm) soft X-ray loop with 
strong emission from two foot points and weak emission 
> 20 keV from the top of the loop. For the June 7 flare, 
the HXR flux > 20 keV (calculated from RHESSI data) 
is £hxr ^9.8x10^'' erg, more than a factor 2000 larger 
than the energy released at high energy. 

SDO/AIA movies at EUV wavelengths show consid- 
erable activity and multiple episodes of loop brighten- 
ing during this flare. SEP proton flux time profiles ob- 
served with GOES above 30, 50 and 100 MeV are plot- 
ted in Figures [l] and [2] There is a significant proton flux 
above 50 MeV starting about 90 min after the March 7 
flare, but no significant flux above 100 MeV is apparent. 
GOES proton fluxes for the June 7 flare increased more 
promptly, were somewhat higher, and extended beyond 
100 MeV. Both flares were associated with fast CMEtE3 
For the March 7 flare, the estimated planc-of-sky veloc- 
ity of the CME is -2000 km s^^, and the CME was pre- 
ceded by an evident shock structurepj while for the June 
7 flare, the CME velocity was measured to be —1000 km 



4. DISCUSSION 

During its first four years, the Fermi-LKT has detected 
7 rays above 100 MeV from more than a dozen solar 
flares, some of which are M-class (see Table II]). This 
suggests that acceleration of electrons and/or protons up 
to several GeV energies may be a common occurrence in 
even modest flares. The highest energy 7 ray recorded 
from the 2011 March 7 flare has an energy of —1 GeV. If 
this is due to pion decay it requires protons of energies 
—5 GeV since the mean energy of 7 rays from i:^ decay 
is typically —1/5 of the proton energy (in the relativistic 
limit). If it is due to electron bremsstrahlung, the bulk of 
the electron distribution responsible for the 7-ray emis- 
sion has an energy of at least 1 GeV, with the exact value 
depending on the electron spectral index. 

Thus, continuous monitoring of the Sun with the 
LAT, in combination with broadband ground-based and 
space-based observations across the electromagnetic and 
charged-particle spectrum, sheds new light on the parti- 
cle acceleration mechanisms in solar flares and perhaps 
elsewhere. 

In this paper we have described the flrst two long- 
duration 7-ray solar flares observed with Fermi-LAT, 
those of 2011 March 7 and June 7. With 13.3 hours 
of 7-ray emission, the March 7 flare was longer than any 
of the flares observed with EGRET. Unlike those earlier 
long-duration 7-ray flares associated only with X-class 



^ ^http://cdaw. gsfc.nasa.gov/CME_list/ 
fiff i,^^'); .1 ^t,-'' ; <,■-- = 



^http: //www. space weather.com/ar chive ■php?vleM=l&day= 
[08&moiith=03&year=201l1 



flares ( Ryan|20"00 ), these two flaress were associated with 
more modest, M-class X-ray flares. The strong detection 
of 7 rays > 100 MeV comprising only a small fraction of 
the energy flux observed in the HXR band demonstrates 
that the sensitivity of the LAT is providing access to a 
new range of solar flare phenomena. Our detailed data 
analysis of these two flares indicates the following: 

1. The > 100 MeV 7 rays clearly originate from the 
Sun and appear to be localized with good confi- 
dence at the active region responsible for the other 
flare activity (seen by RHESSI and SDO). 

2. Although the Sun was not in the Fermi-hKY FOV 
during the impulsive phases of these flares, sim- 
ple extrapolation of the observed light curve of the 
March 7 flare back to the impulsive phase sug- 
gests that significant 7-ray emission may have been 
present then. The comparison between the number 
of detections of impulsive phases ( "Type I" in Ta- 
ble [1]) , and the number of flares with no impulsive 
phase detected ( "Type D" ) , also suggests that this 
may indeed be the case. 

3. The shock front of a CME is known to accelerate 
SEPs, but its contribution to the acceleration of 
the particles that produce 7 rays remains unclear. 
However, as the shock front progresses in interplan- 
etary space, the 7-ray emission cannot occur at the 
shock front itself because densities (<ClO^'^ cm"-^) 
are too low. The accelerated protons and elec- 
trons must be transported to higher densities be- 
low the corona for efficient production of radiation. 
Therefore, the 13.3 hour duration of the gamma- 
ray emission from the 2011 March 7 flare is chal- 
lenging for a shock front acceleration scenario, as 
the CME shock front would have traveled to more 
than about 100 solar radii in this time, making 
transport back to the Sun problematic. Thus, ac- 
celeration in the corona may be more attractive 
as an explanation, especially considering the rela- 
tively large flaring loop seen by RHESSI and the 
extended and complex activity seen by SDO. In this 
case stochastic acceleratio n by turbulence may b e 
the dominant mechanism (Petrosian & Liu 2004). 
The required turbulence can be produced by recon- 
ncction of the magnetic fields in the current sheet 
behind the CME. 

4. The spectral analysis shows that the 7-ray spectra 
cannot be well fit by a simple power law, while 
including an exponential cut-off provides an ac- 
ceptable fit. A model based on decay of pions 
produced by interactions with background parti- 
cles of accelerated protons (and a particles) having 
a power-law spectrum provides an equally accept- 
able fit. We cannot distinguish spectroscopically 
between a model based on bremsstrahlung emis- 
sion from a population of electrons with an intrin- 
sic cut-off and a pion decay model; however, for 
emission mechanisms in the corona-chromosphere, 
electron bremsstrahlung seems unlikely. For typi- 
cal flare densities and magnetic flelds, most (>90%) 
of the > 100 MeV electron energy would appear 
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as sub-mm and far infrared emission via the syn- 
chrotron process and ~50 keV HXRs via inverse- 
Compton scattering of solar optical photons. This 
fraction would be even higher at the lower densities 
expected for the large loop seen by RHESSI. In ad- 
dition, acceleration of electrons to the GeV range 
would require acceleration time scales shorter than 
few seconds to overcome the rapid synchrotron loss 
rate. Whether the leptonic origin of the 7 rays can 
be definitely ruled out by existing observations re- 
quires a more detailed analysis beyond the scope 
of this paper. The hadronic scenario seems more 
plausible and and requires on ly a moderate energy 
input. As mentioned in Sec. 3.4 the total energy 
in long-duration 7 rays above 100 MeV is several 
hundred times less than that of the HXRs observed 
with RHESSI during the short impulsive phase. 
Considering the respective yields of 7 rays by pro- 
tons and HXRs by electrons, we estimate the ratio 
of accelerated proton to electron energy fluxes to 
be < 10~^. So for these M-class flares, the amount 
of energy required to accelerate the protons and 
explain the observations is modest. 



It is common to describe long duration events in 
the framework of the so-called trap-precipitation 
model as was done in the analysis of th e X12.0 
1991 June 11 flare ( [Kanbach et al.||1993| ). If the 
trapping is due to magnetic held convergence and 
if Coulomb collisions are responsible for scattering 
of the particles into the loss cone and their pre- 
cipitation, the trapping efficiency would be higher 
for higher-energy particles. This would result in 
a gradual hardening of the particle spectrum. On 
the contrary, as shown in Table [2] for the pion-decay 



model, we observe a gradual softening of the proton 
spectrum during the 2011 March 7 flare that could 
imply a continuous acceleration process lasting for 
the duration of the emission. 

In summary the Fermi-hAT observations of the loca- 
tions, spectra, and evolution of solar flares in the >100 
MeV energy range have raised interesting issues regard- 
ing acceleration, transport and radiation of particles in 
solar flares. The Fermi LAT Collaboration will continue 
to monitor the Sun through the peak of Solar Cycle 24. 
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